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Abstract 

The contributions of supersymmetric particles in the isospin symmetry violation in B — > K*^ decay 
mode are investigated. The model parameters are adopted from minimal Supergravity with mini- 
mal flavor violation. A complete scan of the mSUGRA parameter space has been performed, using 
the next to leading supersymmetric contributions to the relevant Wilson coefficients. The results 
are compared to recent experimental data in order to obtain constraints on the parameter space. 
We point out that isospin asymmetry can prove to be an interesting observable and imposes severe 
restrictions on the allowed parameter space, in particular for large values of tan [3. The constraints 
obtained with isospin asymmetry also appear as more restricting than the ones from the branching 
ratio of B — > X s j. 

PACS numbers: 11.30.Pb, 12.15.Mm, 12.60.Jv, 13.20.He 



1 Introduction 

The Standard Model (SM) has been very successful to describe the experimental data from accelerator 
physics so far. With new colliders, like the Large Hadron Collider (LHC) or later on the International 
Linear Collider (LLC), becoming operational, the hope is to detect signals which could reveal physics 
beyond the SM, that in turn provide answers to the many theoretical questions left unanswered by 
the SM. 



One of the most motivated scenarios for new physics is generally considered to be Supersymme- 
try (SUSY). In the minimal supersymmetric extension of the Standard Model (MSSM), the large 
number of free parameters makes the phenomenological studies rather complicated. Many studies are 
therefore based on the constrained minimal supersymmetric Standard Model (CMSSM) - often called 
mSUGRA - with the number of parameters reduced to five, corresponding to tjiq (common mass of 
scalar particles at the supersymmetric grand unification scale), m 1 / 2 (universal gaugino mass), A 
(universal trilinear SUSY breaking parameter), together with the sign of the Higgs mixing parameter 
ji and the ratio of the two Higgs vacuum expectation values tan (5. 
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Many studies have been performed to constrain the supersymmetric parameter space, and in par- 
ticular direct and indirect searches for new particles have provided lower bounds on their masses pQ. 
Other constraints come from the cosmological observations of the large scale structures and the cosmic 
microwave background [2], the measurement of the anomalous magnetic moment of the muon {g^ — 2) 
[3], and the study of radiative B meson decays. 

The precision measurements of the radiative B meson decays, which have become possible with the op- 
eration of the B factories and other B-dedicated experiments, have provided exciting opportunities for 
mapping possible routes beyond the SM. One such rare decay mode is the exclusive process B — > K*"f 
and its associated inclusive transition b — > 57, which have been extensively used to constrain new 
physics |H Consequently, a thorough investigation of the branching ratio of these decay modes has 
been instrumental in constraining the parameter space of various models. In this paper, we focus on 
another observable, the isospin asymmetry, and we will show that this observable may even lead to 
more stringent constraints than the branching ratios. 

The isospin asymmetry for the exclusive process B — > K*^ is defined as: 

_ T(B°^K* 7 )-T(B-^K*- 1 ) 
T(B° -» K*°j) + T(B- -> K*~i) ' 1 ' 

with Ao+ obtained from eq.([T]) by using the charge conjugate modes. The most recent data for 
exclusive decays from Belle [6] and Babar [7] point to isospin asymmetries of at most a few percent, 
consistent with zero within the experimental errors: 

A _ = +0.050 ± 0.045(stat.) ± 0.028(syst.) ± 0.024( J R +/0 ) (Babar), (2) 
A 0+ = +0.012 ±0.044(stat.)±0.026(syst.) (Belle), (3) 

where the last error in eq.([2]) is due to the uncertainty in the ratio of the branching fractions of the 
neutral and charged B meson production in Y(4«S) decays. Within the SM, this asymmetry, which is 
due to the non-spectator contributions, has been estimated in the literature, using the QCD factoriza- 
tion approach in Refs. [8] and [9], Brodsky-Lepage formalism [10] . and the perturbative QCD method 
[llj. On the other hand, Ref. [12] deals with the effects of an additional generation of vector-like 
quarks on the isospin symmetry breaking in B — * K*j. We consider here the QCD factorization 
method for our calculations. 



In the following, we first present the general framework of our investigation, followed by an anal- 
ysis of the constraints of the isospin asymmetry on the mSUGRA parameter space in the context 
of minimal flavor violation. In order to compare the constraints due to isospin breaking with those 
obtained from the branching ratios, the results for inclusive branching ratio are reproduced as well. 



2 General framework 

The effective Hamiltonian for b — > s~f transitions reads: 



p=u,c 



d( M ) + c 2 ( M ) of (m) + £ 300 Oi(ji) 



(4) 



i=3 
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where Gp is the Fermi coupling constant, V{j are elements of the CKM matrix, Oj(/i) are the operators 
relevant to B — * K*j and Cj(/i) are the corresponding Wilson coefficients evaluated at the scale \x. 
Since the combination V* s V u b is an order of magnitude smaller than V* s V c b, we can safely neglect the 
u-quark terms. The operators 0« can be listed as follows: 

0{ = s a ^P L ppp^^P L b a , 

Of = SaJ^PLPapp-f^PLbp , 

3 = Sa^Pibg ^ qpl^PLqp , 

9 

4 = Sa^Pzbp qpltiPLQa , (5) 

9 

5 = Sa^PLba^qp^Pnqp , 

9 

06 = Sa^Pzbp ^ qpJij.PRqa , 



7 = -^m b s a a^ v P R b a F^ , 



Os = ^rn b s a a^P R T^b p G% , 

where Pl(Pr) = are projection operators. The electroweak penguin operators are omit- 
ted from the above list as their contributions to the process at hand are negligibly small compared to 
the others. 

The presence of SUSY particles does not introduce new operators in the list, however, the Wilson 
coefficients Cj receive additional contributions from virtual sparticles. 

We use the expressions of the Wilson coefficients at the next to leading order (NLO) in the strong 
coupling constant a s . They are first calculated at the scale fiw = 0(Mw)- The contributions from the 
W boson (SM), the charged Higgs (H) and the charginos (x), as well as the leading tan/3 corrections 
to the W boson and the charged Higgs are considered: 

dfjiw) = C t SM ( m ) + 5C»{n w ) + 8C?(p w ) + 6cf M ' taaf3 \vw) + 5c\ H ' tan/3) ( m ) • (6) 

The contributions from the gluino and neutralino are neglected in our work, as they are known to be 
negligible in the minimal flavor violating scenario [13]. The reason is that within mSUGRA frame- 
work, there exists a strong correlation among various parameters. In particular, the down squarks, 
which appear in the gluino and neutralino loops are much heavier than the stops, the virtual partner 
in the chargino loop. In fact, chargino can be relatively light and since at least one of the stops can 
also be light, their loop results in a considerable contribution. 

The details of the calculation of the Wilson coefficients at the scale fiw are given in Appendix A. 
The Wilson coefficients are then evolved through the renormalization group equations to the scale 
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Hb = 0(mb), at which they can be used to calculate the isospin asymmetry. Further details are given 
in Appendix B. 

Following the method of Ref. [8], one can write the nonspectator isospin symmetry breaking con- 
tribution as A q = b q Ai ea( [, where q is the flavor of the light anti-quark in the B meson and Ai ead 
is the leading isospin symmetry conserving spectator amplitude. To leading order in a s , the main 
contribution to B — > K*^ is from the electromagnetic penguin operator O7: 

Gf 

A lead = -i-^V cb V* s a c 7 {K* 1 \0 7 \B) . (7) 

The factorizable amplitude A[ ea( i is proportional to the form factor T^^ K * which parameterizes the 
hadronic matrix element of O7 to the leading order in Aqcd / m b- The coefficient a^, which is dominated 
by C7, is defined in Appendix C. b q depends on the flavor of the spectator and, in fact, the above 
parameterization leads to a simple expression for the isospin asymmetry in terms of this parameter: 

A _ = Re(b d - b u ) . (8) 

The expression for b q , which is derived in [8] within the QCD factorization method, can be found in 
Appendix C. 

In order to generate the SUSY mass spectrum, as well as the couplings and the mixing matrices, 
we use the Monte Carlo event generator ISAJET-7.74 p3]. We perform scans in the mSUGRA pa- 
rameter space (mo, m^, Aq, sign(/i), tan/3). For any mSUGRA parameter space point, we then 
calculate the isospin asymmetry using eq.©, and compare it to the combined experimental limits of 
eq.© and eq.(|3]). After including the theoretical errors due to the scales and model parameters, we 
allow mSUGRA parameter space points which stand in the 95% confidence level range 

- 0.047 < A _ < 0.093 . (9) 

For comparison, we also perform the calculation of the inclusive branching ratio of B — ► X s ^f following 
Ref. [15] . and allow the mSUGRA parameter space points to be in the 95% confidence level range [lj 

2.33 x 10~ 4 < B(B X sl ) < 4.15 x 10~ 4 . (10) 

The points which result in too small light Higgs masses (i.e. such as m^o < 111 GeV) or which do 
not satisfy the constraints presented in Tabled] are also excluded. 
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Table 1: Lower bounds on sparticle masses in GeV, obtained from pQ. 



Finally, we also examine whether the lightest supersymmetric particle (LSP) is charged. Indeed, the 
LSP is stable when R-parity is conserved, and to be accounted for dark matter, it has to be neutral. 
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On the other hand, if R-parity is violated, then the LSP is not stable and as such cannot be a can- 
didate for the dark matter. In this case, it is possible to have charged LSP with no constraint from 
cosmology. In our results, we have identified the parameter space regions where the LSP is charged 
to indicate the cosmologically disfavored mSUGRA parameters if R-parity is conserved. 

Our results for B — * X s ^f inclusive branching ratio are consistent with those from Ref. [T] and 
from the MicrOMEGAs code [IB] . 

An analysis of our results is presented in the following section. 

3 Constraints from isospin asymmetry 

We perform scans of the mSUGRA parameter space such that mo 6 [0,1000], m\ji G [0,1000], 
tan/3 £ [0,50], Aq G [—1000,1000] and for both signs of For \i < 0, the supersymmetric contri- 
butions to the Wilson coefficients have the same sign as in SM. In this case, the latest experimental 
results are not sufficient to provide constraints on the mSUGRA parameter space. Moreover, fi < 
is disfavored by the (<7« — 2) measurements. For [i > 0, the supersymmetric contributions to the 
Wilson coefficients can have a flipped-sign in comparison to the SM results, leading to a larger isospin 
breaking, and consequently, the experimental data can impose stringent constraints on the mSUGRA 
parameter space. Therefore, in the following we will only present results with /j, > 0. 

An investigation of the (jn-i/2, mo) plane for Aq = is presented in Figure [TJ In this figure, the area 
marked "Isospin" corresponds to the region excluded by the isospin breaking constraints, whereas the 
area marked "B — > X s j n corresponds to the region excluded by the inclusive branching ratio con- 
straints. The "Excluded" area corresponds to the case where at least one of the particle masses does 
not satisfy the constraints of Table [TJ And finally, "Charged LSP" is the cosmologically disfavored 
region when R-parity is conserved. The various colors represent the changing magnitude of the isospin 
asymmetry. 

First, we note that the isospin breaking for a set value of (mi/2,mo) increases with tan/3. Moreover, 
for a fixed tan/3, the asymmetry decreases with larger uiq and m^ 2 . Second, it should be pointed out 
that the constraints from isospin asymmetry are more stringent than the ones from inclusive branch- 
ing ratio. However, for low tan/3 (like tan/3 = 10), constraints from both isospin asymmetry and 
branching ratio are not as restrictive as they are for larger values of tan (5. 

Figure [2] corresponds to the (mirt, mo) plane for Aq = —rriQ. The conventions are the same as in 
the previous figure. Again, we note that the isospin asymmetry is more sensitive to the model param- 
eters than the inclusive branching ratio. A comparison between Figs. [TJ and [2] reveals that the isospin 
symmetry breaking is enhanced by a negative value of Aq so that even for tan/3 as low as 10 it can 
produce appreciable constraints. Nonetheless, the global shapes remain similar. 

The effect of tan (3 on the isospin asymmetry is illustrated in Figure [3j Indeed, the supersymmet- 
ric loop corrections which are proportional to the gluino mass and tan /3 can be quite large at high 
tan/3 limit. This arises from the Hall-Rattazzi-Sarid effect [T7j, and also from the top-quark Yukawa 
coupling [18] . The enhancement of the isospin breaking by tan/3, particularly for smaller values of 
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Figure 1: Constraints on the mSUGRA parameter plane (m^2, mo) f° r = and for different values 
of tan (3. The "Excluded" region in gray corresponds to the sparticle or Higgs search constraints of Ta- 
ble [TJ The light green "Charged LSP" corresponds to a cosmologically disfavored region. The magenta 
U B — > X s ^f" contour delimits the region excluded by the inclusive branching ratio in accordance with 
eq. (llO|) . whereas the yellow "Isospin" region corresponds to the isospin symmetry breaking constraints 
from eq.©. Note that the color scale is different for the first graph with tan/3 = 10. 



m 1 / 2 , is clearly depicted in these graphs. These plots illustrate how stringent the isospin asymmetry 
bounds are at high tan /?, and also reveal the boost of the isospin breaking by a negative value of Aq. 
In fact, the same trend is reported in B s — > decay mode where the branching ratio can improve 

by as much as two orders of magnitude for large values of tan (5 pjJJ [20] . 
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Figure 2: Constraints on the mSUGRA parameter plane (m^, mo) for Aq = —rriQ. The conventions 
for the different regions are the same as in Fig. [TJ 



Figure H] illustrates the sensitivity of the isospin violation to tan/? (and m.1/2) from a different per- 
spective (and somewhat more clear as they are explicit non-contour plots) where the plots are done 
with two Aq values: and — mo- 

To conclude this section, we have shown that the isospin asymmetry can provide stringent constraints 
on the mSUGRA parameter space, and appears to be even more contraining than the inclusive branch- 
ing ratio. 
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tan p tan (3 

Figure 3: Constraints on the mSUGRA parameter plane (tan/3, m^) for mo = 500, with = and 
Aq = —mo. The definitions of the different regions are in the text. 



4 Summary 

In this article, we investigated the possibility for the isospin asymmetry in B — > K*j decay mode to 
be an interesting observable to derive constraints on the supersymmetric parameter space. To obtain 
our results, we calculated the NLO supersymmetric contributions to the isospin asymmetry, using the 
effective Hamiltonian approach within the QCD factorization method and considering the minimal 
flavor violation. The mSUGRA parameter space was scanned, and the resulting isospin asymmetry 
for each point was compared to the experimental data from Babar and Belle. Our main conclusion of 
this comparison is that, provided /i > 0, the isospin asymmetry appears to be a powerful observable to 
constrain the mSUGRA parameter space producing even more stringent restrictions than the inclusive 
branching ratio. Among the different parameters, the values of and tan (3 appear to be restricted 
more significantly by the isospin symmetry breaking constraints. 

In this work we considered the mSUGRA model which has the advantage of having a fewer number 
of free parameters. However, mSUGRA's assumptions are in fact very strong and therefore extending 
this study to other SUSY models can be of interest. 

To conclude, with more accurate experimental data, we can hope the isospin asymmetry could re- 
veal to be a very valuable observable to explore the supersymmetric parameter space. 
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Figure 4: Isospin asymmetry versus mi/2 and tan/3 for ^0 = and Aq 
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Appendix A: Wilson coefficients at /i w 

The Wilson coefficients follow a perturbative expansion: 



4tt 



where the a s evolution writes pQ: 

-in 



a s (p) 



with 



/3 ln(M 2 /A(»/) 



ft ln[ln( M 2 /A("/) 2 )] 



PI 



lm/^/A ; .«iuV/A •■*) i v ln[ln(/i2/A( " /)2)] 2 



2 38 , „ 5033 325 „ 

/3 = H--n / , ft = 102-— and = 2857 - -^n f + — n) 



(Al) 



&A) 5 
2/3 2 4 



(A2) 
(A3) 



3 ' 3 J ' ^ 9 J 27 

being the number of flavors. A( n /) is a dimensional parameter depending on the number of flavors. 
The numerical values of A's in table [3] are based on the input a s (Mz) = 0.1172. 

The main contributions to the Wilson coefficients are classified into three categories: 1) those from 
the Standard Model, 2) charged Higgs contributions and 3) chargino contributions. The details of 
each contributing term are given in the following. 



A.l Standard Model contributions 

The Standard Model contributions to the Wilson coefficients are adopted from Ref. [21] . At leading 
order (LO), they read: 



C ^°\ m )=cr w ^w) = { 



SM(0) , 



for 1 = 1,3,4,5,6 

1 for i = 2 



where 



xtw 



F^(x tw ) for i = 7,8 , 



F 7 (1) (x) 



F 8 (1) (x) 



x(7-5x-8x 2 ) x 2 (3x-2), 
24(x-l) 3 4(x-l) 4 ' 

x(2 + 5x — x 2 ) 3x 2 



8(x-l) 3 4(x-l) 4 
The NLO top quark running mass at a scale fi is given by [TJ [52] : 



lnx . 



m t (p) = m t {m t ) 





10 
2/3 


_a s (m t )_ 





a s {m t ) 7o m HJ_ _ M ( <*.(/*) _ 
4vr 2f3 V7o m M \« S K) 



(A4) 



(A5) 



(A6) 



(A7) 
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and 



m t (m t ) = m t 



_ 4a s (m t ) 



3 vr 

Here m t is the pole mass of the top quark. (3q and j3\ are defined in eq. (|A3P and: 



76" 



The NLO corrections can be written as 1211: 



404 40 



n SM{l), \ . 



7i 



15 + 61n 




-n f 



(A8) 



(A9) 



for 



t = 1 



for i = 2, 3, 5, 6 



2 2 , & 



for 



t = 4 



(A10) 



3 ' 3" M*y 
Gi(x t w) + Ai(x t w)fo-T7T~ for i = 7,8 , 



Mi 



where 
E(x) 

G 7 (x) 



w 



x(-18 + llx + a; 2 ) x 2 (15 - 16x + 4x 2 ) . 2, 

; r5 1 ; m In x In x , 
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(All) 
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9(x - If 



Li 2 1 



81 (s - l) 5 



lnx + 



-28x 2 + 46x 3 + 6x 4 
3(x - l) 5 



ln 2 x , (A12) 



Gh(x) 



A 7 (x) 



A s fx) 



-508 + 610x- 28209x 2 - 14102x 3 + 737x 4 x + 41x 2 + 40x 3 - 4x 4 L _ / 1 
1296(x - l) 4 h 6(x - l) 4 12 V ~ x 
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H ; rs In x H ; r= — in x , 
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140 - 902x - 1509x 2 - 398x 3 + 77x 4 x 2 (31 + 17 x) 



lnx 
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+ 



2(x - 1) 



and where Li2 is the usual dilogarithm function Li2(x) 



lnx 



ln(l - t) 



(A13) 
(A14) 
(A15) 



dt 

o t 



A. 2 Charged Higgs contributions 

At the LO, the relevant charged Higgs contributions to the Wilson coefficients are given by [21]: 



l -F% (x tH ± ) - A u A d F% (x tH ± ) , 



(A16) 
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with 



where 



A __J_ 



1 , m 2 t {^ w ) 

and x tH ± — 



F 7 (2) (x) 
if } (x) 



Arf tan (5 

x(3 — 5x) x(3x — 2) 



(A17) 
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12(x-l) 2 6(x-l) 3 

x(3 — x) x 
4(x - l) 2 ~ 2(x - 1) E 



lnx , 



lnx 



(A18) 



and F 7 (1 8 } are defined in eq. (|A6p . 

At the NLO, the charged Higgs contributions are: 



Sdp{ji w ) = G«{x tH± ) + ^(x tH± )\n^--E H (x tH± ) 
8C$\nw) = Gf(x tH± ) + A^(x tH± )\n^-^E H (x tH± ) 



(A19) 
(A20) 



with 
Gf(x) 



4 



4(-3 + 7x- 2x 2 



3(x - l) s 



Li 2 1 + 
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1 
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9(x - l) 5 

21 - 47x + 8x 2 2(-8 + 14x + 3x 2 ) 



108(x - l) 4 



(A21) 



+ 



lnx 



(x - l) 3 ! (x - l) 4 

-31 - 18x + 135x 2 - 14x 3 x(14 - 23x - 3x 2 ) 



6(x- l) 4 
-36 + 25x - 17x 2 
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(x - iy 



lnx 



(A22) 



2(x - l) 3 

3(143 - 44x + 29x 2 
8(x - l) 3 
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A. 3 Chargino contributions 

2 

Tfl • 

In the following, we use the notation Xj,- = — % . The masses of the sparticles are assumed such that 
for the squarks mg, < rrig , for the charginos m + < m+, and for the neutralinos m v o < m v o < 

^ y Xl X2 Xl A2 

m^o < m^o . The relevant chargino contributions read |23[ [23] : 

- -EE {|l^#F«(x faf ) + r;,r^FS»(: Cfix± )) 

fc=l i=l I *fc Xj J 

+ E(^P»l#^(^f)+f;A^F«( w )} , (A26) 

i=l I 912 Xi ) 

where fi s is the SUSY scale, and mq 12 is the common mass of the up and charm squarks, which we 
consider identical (mg 12 ~ ~ mg). Moreover, we have 



r'ii = -s J, 1 *. 3 .. ^ > (A27) 



\/2cos/3(l + e£tan/3) ' 



where £7 and V are the chargino mixing matrices, following the diagonalizing convention: 

^Uv^v^cos/? ^ J F ' (A28) 
and is the squark g mixing matrix such as: 

Da=(™ 9 J , (A29) 



and eb, which will be given below, is a two loop SUSY correction, whose effects are enhanced by factors 
of tan (5. 

Ty and T'- are obtained from Tij and V- by replacing the matrix by the unity matrix. The 
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functions ivg(a;) are given by [23] : 



,(3)/ n _ (5 -7a) x(3x-2) 
7 W 6(x-l) 2 3(x-l) 3 



The value of the chargino contributions at the scale \xy/ is computed using: 



= ^ 3/30 6C*(ii s ) + - U 3/3(1 - r? s 3 "° ] 5C 8 % S ) , (A31) 



14 
M 7 " 



= n s * P °5C*(n s ) , (A32) 
where rj s = a s (n s )/a s (fiw) and f3' Q = —7, which corresponds to six active flavors. 
In the following, we adopt the notations: 

COS#g = Ajll = Cq , SUlBq = Dq 2 l = Sq . (A33) 

The leading tan/3 corrections are contained in the following formulas for e&, e' b and e' t , which are 
evaluated at scale fi s [TBJ E3] : 



P6 



2a s {n s ) A b /t<m(3 - \i . Vt(Ps) \^ tt ^/tan(3-A t 



-^Ks'^ + T^r 2^ ^ — ^~ H Hxf' x hxf^ 



3vr m a 619 629 16 tt 2 ^ m ± 

9 i=l,2 x< 



q(M z )/iM 2 
4 sin 2 6>h/tt 

-? / n/f2 ,,2 \ «2 / A/f 2 ,,2 \ r- 2 / ,,2 \ „? / A#2 ,,2 



+^4^= x (A34) 



M ±iU-fLff M JfL] +^Lh Mi j?L\+JLtt Mi J£ 



2 I 2' 2 I 2 " I 2' 2l'n2 il l 2' 2 I 1 n 2 Ji I 2' 2 

m~ \ m~ ] m~ \ m~ m~ } 2m- \ m- m- / 2m- \ m- m- 

tl \ tl tl / 12 V *2 t2 / 61 \ 61 61 / 62 \ 02 02 . 



where A g is the trilinear coupling of the quark q. y q and y q are the ordinary and supersymmetric 
Yukawa couplings of the quark q respectively. The function H is defined by: 

H{x 1 y) = T - + - r . (A35) 

(l-x)(x-y) (l-y)(y-x) 

Please note that we neglect the neutralino mixing matrices and we assume that the chargino masses 
are given by fi and Mi. 
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, _ 2a s (fx s ) A b / tan 0-fj, 
e b [t) — — : x 



3?r 



44 H ( X ti~9> X b 2 g) + 44 H ( X hg> X ~hg) + 44 H ( X i2~g> X b 2 g) + 44 H ( X h~g> X ~hg) 

A t — \xj tan (3 



16 vH 



ii' 



i=l 



c t c fe ff ( x tox?> ^x?) + c t - Sg H(x hx o, x~ b2X o) + s~ t c h H(x kx o,x~ bix o) + s f s~ b H(x~ tlx o, x-^ 



| a{Mz)nM 2 ^ 
4 sin 2 # W 7T 



M| M 2 



s 2 



2 ' 2 

> 01 01 . 



+ 



-H 



Ml fj? 



(A36) 

N i3 



m 



In 



2 ' 2 
TO- m~ 
02 02 . 



ti 



to- ' to- 



+ 



2m? 

*2 



if 



M| ^ 2 



2 ' 2 
*2 12 , 



iV, in the above formula, represents the neutralino mixing matrix. The last correction reads: 



m-n 



+ 



2a s fi + At I tan (3 
1/6 (M* 



4 H ( x hg> x ~sg) + 4 H ( x hgi x ~s~9) 



(A37) 



16 vH 



i=l 



/// tan/? 
m v o 



C H y? > x 6, v« ) + c ? 5 S ff ( x h y9 ' ^ v° ) 



2„2 



iV i3 



The SM and charged Higgs contributions at the nw scale are affected by e , e( and e' t as the following: 



(SM.tan/J) 
7,8 



11" / 



fa> ~ tan /? p (2)/ v 
l + €fe tan/? F ^ Xtw) ' 



(A38) 
(A39) 



1 + e b tan /3 

Finally, the complete Wilson coefficients Ci°^ are found by adding the different contributions: 

cftW) = Crj^W) + ^rTW) + fCftfar) + 5Cg A/ ' tan/?) ( m ) + ^J tan/3 V^) (A40) 
cftW) = C 7 S f (1) ( MW ) + <5Cf 8 (1) ( m ) • (A41) 



Appendix B: Wilson coefficients at //& 

The Wilson coefficients at the lower scale = 0{m b ) can be written as [221 125j : 

CM) = Cf\» b ) + ^^0*) + • • • , (Bl) 
15 



where, for j = 1 • • • 6: 

i=3 i=3 



with 



8 

cf (/x 6 ) = ^cf\ m ) + C^\ m )Y,hri ai ■ 



i=l 



The next to leading coefficient is given by 



( 297664 is 7164416 w 256868 37 6698884 39 \ t 0)/ 

+ 7723 T7 23 -I 7723 T7 23 Co [Uw) 

V 14283 357075 ' 14283 ' 357075 ' / 8 VP ; 

_ 8 

+ W ( ?7§ " C 7° ] ^ + YtiwE&w) +fi + 9iV)v a * + AC«( W 



i=l 

,2 



f 2 \ 

where in the M5" scheme AC^/ij) = ( -ej + 6/j ) rf i+l In . The numbers aj, fc^, e 

i= i V 3 / M w 

§ij , hi, hi, ei, gi and k are gathered in Table [2]0. 



Appendix C: Isospin asymmetry 

To leading order the isospin asymmetry Ao- is given by [8]: 

A _ = Re(b d - b u ) . 

The spectator-dependent coefficients b q can be written as: 

12tt 2 / b Q q f f£, f K *m K . \ 

m b Tf n aj \ m b 6\ B m B J 

1 As we follow the operator basis from |25j . we also use the numerical values of this article. 
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i 


1 


2 


3 


4 


5 


6 


7 


8 




14 

23 


16 

23 


6 

23 


12 

23 


0.4086 


-0.4230 


-0.8994 


0.1456 


hi 








1 

2 


1 
2 














en 


























fli 








0.8136 


0.7142 














9ii 








1.0197 


2.9524 














hi 








i 

2 


i 

2 








































hi 








0.8136 


-0.7142 














92% 








1.0197 


-2.9524 














hi 








l 

14 


l 

6 


0.0510 


-0.1403 


-0.0113 


0.0054 


e3i 














0.1494 


-0.3726 


0.0738 


-0.0173 


hi 








-0.0766 


-0.1455 


-0.8848 


0.4137 


-0.0114 


0.1722 


93i 








-0.1457 


-0.9841 


0.2303 


1.4672 


0.0971 


-0.0213 


k^i 








i 

14 


i 

6 


0.0984 


0.1214 


0.0156 


0.0026 
















0.2885 


0.3224 


-0.1025 


-0.0084 


hi 








-0.2353 


-0.0397 


0.4920 


-0.2758 


0.0019 


-0.1449 


94i 








-0.1457 


0.9841 


0.4447 


-1.2696 


-0.1349 


-0.0104 


hi 














-0.0397 


0.0117 


-0.0025 


0.0304 
















-0.1163 


0.0310 


0.0162 


-0.0975 


hi 








0.0397 


0.0926 


0.7342 


-0.1262 


-0.1209 


-0.1085 


95i 














-0.1792 


-0.1221 


0.0213 


-0.1197 


hi 














0.0335 


0.0239 


-0.0462 


-0.0112 
















0.0982 


0.0634 


0.3026 


0.0358 


hi 








-0.1191 


-0.2778 


-0.5544 


0.1915 


-0.2744 


0.3568 


9m 














0.1513 


-0.2497 


0.3983 


0.0440 


hi 


2.2996 


-1.0880 


3 

7 


l 

14 


-0.6494 


-0.0380 


-0.0185 


-0.0057 


hi 


0.8623 











-0.9135 


0.0873 


-0.0571 


0.0209 


e>i 


4661194 

816831 


8516 
2217 








-1.9043 


-0.1008 


0.1216 


0.0183 


fi 


-17.3023 


8.5027 


4.5508 


0.7519 


2.0040 


0.7476 


-0.5385 


0.0914 


9i 


14.8088 


-10.8090 


-0.8740 


0.4218 


-2.9347 


0.3971 


0.1600 


0.0225 


h 


0.5784 


-0.3921 


-0.1429 


0.0476 


-0.1275 


0.0317 


0.0078 


-0.0031 



Table 2: Useful numbers. 



The dimensionless functions K\ and K2 q are given by: 
K\ = - 



+■ 



N 4vr 



m B J 



3 fi b 3 



- G±(x cb ) 



(C3) 



+ ri 
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A" 



2q 



VusVub 



v: s v cb 

Cf cx s ((j, b ) 



C 2 (fx b ) + 



N 



Cz{nb 
N 



(C4) 



+ 



N 4?r 

where N = 3 and Cf = 4/3 are color factors, and 



(Mm,) ( y| In ^ • ^ // . (>,.,,) ) + r 2 



r-2 



3 G^ft) + 3 n/(C 4 (M6) + C 6 (m 6 )) - 8 (JVC 6 (/i ft ) + C 5 (/4»)) 



A ± ln^ + ... 
MO 



44 4 

-g C^Mfe) - 3 n/(C 4 (/i 6 ) + C 6 (Mb)) 



ln^ + 

Mo 



(C5) 



here nj = 5, and = 0(rrib) is an arbitrary normalization scale. r*i and r2 are neglected in our 
calculations. The pole mass of the quarks can be deduced from the running quark mass at fh q pQ: 



m„ 



m q (m q ) < 1 + 



4a s (m g ) 
3vr 



+ 



1 _ , _ 9 , ) + 13.4434 

\ 3 m q (m q ) 



a s (m q ) 



TT 



(C6) 







"a s (m 9 )~ 




7T 





+ [0.6527nJ - 26.655n /; + 190.595] 

where the sum over k extends over the nj t flavors of the quarks lighter than the quark q. The func- 
tions F±, G±(x c b), H±(x c b) and X±_ are convolution integrals of hard-scattering kernels with the meson 
distribution amplitudes, their values are given in Table El The parameter X = ln(ms/A/j) (1 + ge llp ) 
in this table parameterizes the logarithmically divergent integral Jq dx/(l — x). We have evaluated 
the theoretical uncertainty by allowing g < 1 and the phase <p to be arbitrary. A/j ~ 0.5 GeV is a 
typical hadronic scale. 



The coefficient a 7 is given by [9p: 



a c 7 (iT 7 ) = C 7 ( M6 )+ 



47T 



C 2 {u. b )G 2 {x cb )+C 8 {u. b )G 8 ) + 



Air 



C 2 (fi h )H 2 (x cb )+C 8 (^ h )H 8 

(C7) 



in which fih = V^-h^b, and 



^ / s 104 . Hb . , v ^ 8 

k 2 (x c fc) = — —In h 32(a?c6j , &8 = o ln 1" 98 

21 mb o nib 



(C8) 



with 



^(x) = -(48 + 30«vr-57r 2 -2i7r 3 -36C3 + (36 + 6i7r-97r 2 )lnx + (3 + 6i7r)ln 2 x + ln ; 



X X 



+ - (18 + 2vr 2 - 2ivr 3 + (12 - 6vr 2 ) ln x + 6in ln 2 x + In 3 a;) x 1 



(C9) 



+ ^=( - 9 + 112nr - 14?r 2 + (182 - 48ivr) lnx - 1261n 2 x) 

— ■ I 

11 2vr 2 2m 



2 ^ 3 833 20ivr 8vr 2 3/2 
" X ~ 162 " ~W + ~9~ X 



+ 



3 ' 



(CIO) 



A comparison has also been performed with the results of 
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CKM and coupling constant parameters 


V us 




\v ub /v cb \ 


Re(V: s V ub /V c * s V cb ) 


A (4) 


A (5) 


A (6) 


0.22 


0.041 ±0.02 


0.085 ± 0.025 


0.011 ± 0.005 


0.277 GeV 


0.200 GeV 


0.085 GeV 



Parameters related to the B mesons 


TUB 


fB 


Xb 


H 2 (x cb ) 




5.28 GeV 


200 ± 20 MeV 


(350 ± 150) MeV 


(-0.27 ± 0.06) + (-0.35 ± 0.10)* 


0.70 ± 0.07 



Parameters related to the K* meson 


m K * 


fx* 


fa 


jiB-^K* 


892 MeV 


226 ± 28 MeV 


175 ± 9 MeV 


0.30 ±0.05 



Parameters related to the convolution integrals 


F± 


G±(xcb) 


H±(x cb ) 


x± 


1.21 ±0.06 


(2.82 ±0.20) + (0.81 ±Q.23)« 


(2.32 ±0.16) + (0.50 ±0.18)i 


(3.44 ± 0.47) X - (3.91 ± 1.08) 



Quark and W-boson masses 


m b {m b ) 


m c (m c ) 


m t 


M w 


4.20 ± 0.07 GeV 


1.25 ± 0.09 GeV 


174.2 ± 3.3 GeV 


80.4 GeV 



Table 3: The numerical values of the used parameters. 



where C3 ~ 1.2020569 and x 



in 



cb 



|- . We have also: 



H 2 {x) 



2vr 2 f B fk* f 1 ^bi(Q f\ , n ^ , v 

j — 2 I at, — - — / dvh{l — v,x)$±(v) 
i r n B Jo ? Jo 



3.Y l( : i- 



where h(u, s) is the hard-scattering function given by: 



u \ Ax 
h(u,x) = — 



u- 



Li2 



V 



u — 4x + is 



+ Li 2 



u 



I 



1 + 



u — 4x + ie 



I 



2 

u 



(Cll) 



(C12) 



<E>x is the light-cone wave function with transverse polarization and $bi is a distribution amplitude 
of the B meson involved in the leading-twist projection. In a first approximation, $x can be reduced 
to its asymptotic limit <&±(x) = 6x(l — x) [27] . Finally, we can write: 



mT B-+K* m 2 B J £ J Q 



(C13) 



The first negative moment of &bi is parameterized by the quantity Xb such as / d£ — = — — 

Jo £ X B 

The values of the different parameters can be found in Tabled 
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Using these relations altogether, it is then possible to calculate the isospin asymmetry from eq, (|Cip . 



For the computation of the inclusive branching ratios, we also used the relations contained in the 
appendixes, together with those of Ref. [T5] . 
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